A second hotspot for pathogenic exon-skipping variants in CDC45.

Biallelic pathogenic variants in CDC45 are associated with Meier-Gorlin syndrome with craniosynostosis (MGORS type 7), which also includes short stature and absent/hypoplastic patellae. Identified variants act through a hypomorphic loss of function mechanism, to reduce CDC45 activity and impact DNA replication initiation. In addition to missense and premature termination variants, several pathogenic synonymous variants have been identified, most of which cause increased exon skipping of exon 4, which encodes an essential part of the RecJ-orthologue's DHH domain. Here we have identified a second cohort of families segregating CDC45 variants, where patients have craniosynostosis and a reduction in height, alongside common facial dysmorphisms, including thin eyebrows, consistent with MGORS7. Skipping of exon 15 is a consequence of two different variants, including a shared synonymous variant that is enriched in individuals of East Asian ancestry, while other variants in trans are predicted to alter key intramolecular interactions in α/ß domain II, or cause retention of an intron within the 3'UTR. Our cohort and functional data confirm exon skipping is a relatively common pathogenic mechanism in CDC45, and highlights the need for alternative splicing events, such as exon skipping, to be especially considered for variants initially predicted to be less likely to cause the phenotype, particularly synonymous variants.

Widespread repression of anti-CRISPR production by anti-CRISPR-associated proteins.

Many bacteria use CRISPR-Cas systems to defend against invasive mobile genetic elements (MGEs). In response, MGEs have developed strategies to resist CRISPR-Cas, including the use of anti-CRISPR (Acr) proteins. Known acr genes may be followed in an operon by a putative regulatory Acr-associated gene (aca), suggesting the importance of regulation. Although ten families of helix-turn-helix (HTH) motif containing Aca proteins have been identified (Aca1-10), only three have been tested and shown to be transcriptional repressors of acr-aca expression. The AcrIIA1 protein (a Cas9 inhibitor) also contains a functionally similar HTH containing repressor domain. Here, we identified and analysed Aca and AcrIIA1 homologs across all bacterial genomes. Using HMM models we found aca-like genes are widely distributed in bacteria, both with and without known acr genes. The putative promoter regions of acr-aca operons were analysed and members of each family of bacterial Aca tested for regulatory function. For each Aca family, we predicted a conserved inverted repeat binding site within a core promoter. Promoters containing these sites directed reporter expression in E. coli and were repressed by the cognate Aca protein. These data demonstrate that acr repression by Aca proteins is widely conserved in nature.

Analysis of Fungal Genomes Reveals Commonalities of Intron Gain or Loss and Functions in Intron-Poor Species.

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron-exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.

Quantitative analysis of the splice variants expressed by the major hepatitis B virus genotypes.

Hepatitis B virus (HBV) is a major human pathogen that causes liver diseases. The main HBV RNAs are unspliced transcripts that encode the key viral proteins. Recent studies have shown that some of the HBV spliced transcript isoforms are predictive of liver cancer, yet the roles of these spliced transcripts remain elusive. Furthermore, there are nine major HBV genotypes common in different regions of the world, these genotypes may express different spliced transcript isoforms. To systematically study the HBV splice variants, we transfected human hepatoma cells, Huh7, with four HBV genotypes (A2, B2, C2 and D3), followed by deep RNA-sequencing. We found that 13-28 % of HBV RNAs were splice variants, which were reproducibly detected across independent biological replicates. These comprised 6 novel and 10 previously identified splice variants. In particular, a novel, singly spliced transcript was detected in genotypes A2 and D3 at high levels. The biological relevance of these splice variants was supported by their identification in HBV-positive liver biopsy and serum samples, and in HBV-infected primary human hepatocytes. Interestingly the levels of HBV splice variants varied across the genotypes, but the spliced pregenomic RNA SP1 and SP9 were the two most abundant splice variants. Counterintuitively, these singly spliced SP1 and SP9 variants had a suboptimal 5' splice site, supporting the idea that splicing of HBV RNAs is tightly controlled by the viral post-transcriptional regulatory RNA element.

Discovery of multiple anti-CRISPRs highlights anti-defense gene clustering in mobile genetic elements.

Many prokaryotes employ CRISPR-Cas systems to combat invading mobile genetic elements (MGEs). In response, some MGEs have developed strategies to bypass immunity, including anti-CRISPR (Acr) proteins; yet the diversity, distribution and spectrum of activity of this immune evasion strategy remain largely unknown. Here, we report the discovery of new Acrs by assaying candidate genes adjacent to a conserved Acr-associated (Aca) gene, aca5, against a panel of six type I systems: I-F (Pseudomonas, Pectobacterium, and Serratia), I-E (Pseudomonas and Serratia), and I-C (Pseudomonas). We uncover 11 type I-F and/or I-E anti-CRISPR genes encoded on chromosomal and extrachromosomal MGEs within Enterobacteriaceae and Pseudomonas, and an additional Aca (aca9). The acr genes not only associate with other acr genes, but also with genes encoding inhibitors of distinct bacterial defense systems. Thus, our findings highlight the potential exploitation of acr loci neighborhoods for the identification of previously undescribed anti-defense systems.

Purple haze: Cryptic purple sequestrate Cortinarius in New Zealand.

CORTINARIUS: is a species-rich ectomycorrhizal genus containing taxa that exhibit agaricoid or sequestrate basidiome morphologies. In New Zealand, one of the most recognizable and common Cortinarius species is the purple sequestrate fungus, C. porphyroideus. We used genome skimming of the almost 100-y-old type specimen from C. porphyroideus to obtain the nuc rDNA internal transcribed spacer region ITS1-5.8S-ITS2 (ITS barcode) and partial nuc rDNA 28S (28S) sequences. The phylogenetic position of C. porphyroideus was established, and we found that it represents a rarely collected species. Purple sequestrate Cortinarius comprise multiple cryptic species in several lineages. We describe four new species of Cortinarius with strong morphological similarity to C. porphyroideus: Cortinarius diaphorus, C. minorisporus, C. purpureocapitatus, and C. violaceocystidiatus. Based on molecular evidence, Thaxterogaster viola is recognized as Cortinarius violaceovolvatus var. viola. These species are associated with Nothofagus (southern beech) and have very similar morphology to C. porphyroideus but are all phylogenetically distinct based on molecular data.

Publisher Correction: Targeting of temperate phages drives loss of type I CRISPR-Cas systems.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A ribosomal sequence-based oil sensitivity index for phytoplankton groups.

Species-level variability has made it difficult to determine the relative sensitivity of phytoplankton to oil and mixtures of oil and dispersant. Here we develop a phytoplankton group sensitivity index using ribosome sequence data that we apply to a mesocosm experiment in which a natural microbial community was exposed to oil and two oil-dispersant mixtures. The relative sensitivity of four phytoplankton taxonomic groups, diatoms, dinoflagellates, green algae, and Chrysophytes, was computed using the log of the ratio of the number of species that increase to the number that decrease in relative abundance in the treatment relative to the control. The index indicates that dinoflagellates are the most sensitive group to oil and oil-dispersant treatments while the Chrysophytes benefit under oil exposure compared to the other groups examined. The phytoplankton group sensitivity index can be generally applied to quantify and rank the relative sensitivity of diverse microbial groups to environmental conditions and pollutants.

Targeting of temperate phages drives loss of type I CRISPR-Cas systems.

On infection of their host, temperate viruses that infect bacteria (bacteriophages; hereafter referred to as phages) enter either a lytic or a lysogenic cycle. The former results in lysis of bacterial cells and phage release (resulting in horizontal transmission), whereas lysogeny is characterized by the integration of the phage into the host genome, and dormancy (resulting in vertical transmission)1. Previous co-culture experiments using bacteria and mutants of temperate phages that are locked in the lytic cycle have shown that CRISPR-Cas systems can efficiently eliminate the invading phages2,3. Here we show that, when challenged with wild-type temperate phages (which can become lysogenic), type I CRISPR-Cas immune systems cannot eliminate the phages from the bacterial population. Furthermore, our data suggest that, in this context, CRISPR-Cas immune systems are maladaptive to the host, owing to the severe immunopathological effects that are brought about by imperfect matching of spacers to the integrated phage sequences (prophages). These fitness costs drive the loss of CRISPR-Cas from bacterial populations, unless the phage carries anti-CRISPR (acr) genes that suppress the immune system of the host. Using bioinformatics, we show that this imperfect targeting is likely to occur frequently in nature. These findings help to explain the patchy distribution of CRISPR-Cas immune systems within and between bacterial species, and highlight the strong selective benefits of phage-encoded acr genes for both the phage and the host under these circumstances.

Genome-wide correlation analysis suggests different roles of CRISPR-Cas systems in the acquisition of antibiotic resistance genes in diverse species.

CRISPR-Cas systems are widespread in bacterial and archaeal genomes, and in their canonical role in phage defence they confer a fitness advantage. However, CRISPR-Cas may also hinder the uptake of potentially beneficial genes. This is particularly true under antibiotic selection, where preventing the uptake of antibiotic resistance genes could be detrimental. Newly discovered features within these evolutionary dynamics are anti-CRISPR genes, which inhibit specific CRISPR-Cas systems. We hypothesized that selection for antibiotic resistance might have resulted in an accumulation of anti-CRISPR genes in genomes that harbour CRISPR-Cas systems and horizontally acquired antibiotic resistance genes. To assess that question, we analysed correlations between the CRISPR-Cas, anti-CRISPR and antibiotic resistance gene content of 104 947 reference genomes, including 5677 different species. In most species, the presence of CRISPR-Cas systems did not correlate with the presence of antibiotic resistance genes. However, in some clinically important species, we observed either a positive or negative correlation of CRISPR-Cas with antibiotic resistance genes. Anti-CRISPR genes were common enough in four species to be analysed. In Pseudomonas aeruginosa, the presence of anti-CRISPRs was associated with antibiotic resistance genes. This analysis indicates that the role of CRISPR-Cas and anti-CRISPRs in the spread of antibiotic resistance is likely to be very different in particular pathogenic species and clinical environments. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Bioinformatic evidence of widespread priming in type I and II CRISPR-Cas systems.

CRISPR-Cas systems provide bacteria and archaea with adaptive immunity against invading genetic elements, such as plasmids, bacteriophages and archaeal viruses. They consist of cas genes and CRISPR loci, which store genetic memories of previously encountered invaders as short sequences termed spacers. Spacers determine the specificity of CRISPR-Cas defence and immunity can be gained or updated by the addition of new spacers into CRISPR loci. There are two main routes to spacer acquisition, which are known as naïve and primed CRISPR adaptation. Naïve CRISPR adaptation involves the de novo formation of immunity, independent of pre-existing spacers. In contrast, primed CRISPR adaptation (priming) uses existing spacers to enhance the acquisition of new spacers. Priming typically results in spacer acquisition from locations near the site of target recognition by the existing (priming) spacer. Primed CRISPR adaptation has been observed in several type I CRISPR-Cas systems and it is potentially widespread. However, experimental evidence is unavailable for some subtypes, and for most systems, priming has only been shown in a small number of hosts. There is also no current evidence of priming by other CRISPR-Cas types. Here, we used a bioinformatic approach to search for evidence of priming in diverse CRISPR-Cas systems. By analysing the clustering of spacers acquired from phages, prophages and archaeal viruses, including strand and directional biases between subsequently acquired spacers, we demonstrate that two patterns of primed CRISPR adaptation dominate in type I systems. In addition, we find evidence of a priming-like pathway in type II CRISPR-Cas systems.

Response of natural phytoplankton communities exposed to crude oil and chemical dispersants during a mesocosm experiment.

During the 2010 Deepwater Horizon oil spill, the chemical dispersant Corexit was applied over vast areas of the Gulf of Mexico. Marine phytoplankton play a key role in aggregate formation through the production of extracellular polymeric materials (EPS), an important step in the biological carbon pump. This study examined the impacts of oil and dispersants on the composition and physiology of natural marine phytoplankton communities from the Gulf of Mexico during a 72-hour mesocosm experiment and consequences to carbon export. The communities were treated using the water accommodated fraction (WAF) of oil, which was produced by adding Macondo surrogate oil to natural seawater and mixed for 24 h in the dark. A chemically enhanced WAF (CEWAF) was made in a similar manner, but using a mixture of oil and the dispersant Corexit in a 20:1 ratio as well as a diluted CEWAF (DCEWAF). Phytoplankton communities exposed to WAF showed no significant changes in PSII quantum yield (Fv/Fm) or electron transfer rates (ETRmax) compared to Control communities. In contrast, both Fv/Fm and ETRmax declined rapidly in communities treated with either CEWAF or DCEWAF. Analysis of other photophysiological parameters showed that photosystem II (PSII) antenna size and PSII connectivity factor were not altered by exposure to DCEWAF, suggesting that processes downstream of PSII were affected. The eukaryote community composition in each experimental tank was characterized at the end of the 72 h exposure time using 18S rRNA sequencing. Diatoms dominated the communities in both the control and WAF treatments (52 and 56% relative abundance respectively), while in CEWAF and DCEWAF treatments were dominated by heterotrophic Euglenozoa (51 and 84% respectively). Diatoms made up the largest relative contribution to the autotrophic eukaryote community in all treatments. EPS concentration was four times higher in CEWAF tanks compared to other treatments. Changes in particle size distributions (a proxy for aggregates) over time indicated that a higher degree of particle aggregation occurred in both the CEWAF and DCEWAF treatments than the WAF or Controls. Our results demonstrate that chemically dispersed oil has more negative impacts on photophysiology, phytoplankton community structure and aggregation dynamics than oil alone, with potential implications for export processes that affect the distribution and turnover of carbon and oil in the water column.

Prediction and diversity of tracrRNAs from type II CRISPR-Cas systems.

Type II CRISPR-Cas9 systems require a small RNA called the trans-activating CRISPR RNA (tracrRNA) in order to function. The prediction of these non-coding RNAs in prokaryotic genomes is challenging because they have dissimilar structures, having short stems (3-6 bp) and non-canonical base-pairs e.g. G-A. Much of the tracrRNA is involved in base-pairing interactions with the CRISPR RNA, or itself, or in RNA-protein interactions with Cas9. Here we develop a new bioinformatic tool to predict tracrRNAs. On an experimentally verified test set the algorithm achieved a high sensitivity and specificity, and a low false discovery rate (FDR) on genome analysis. Analysis of representative RefSeq genomes (5462) detected 275 tracrRNAs from 165 genera. These tracrRNAs could be grouped into 15 clusters which were used to build covariance models. These clusters included Streptococci and Staphylococci tracrRNAs from the CRISPR-Cas9 systems which are currently used for gene editing. Compensating base changes observed in the models were consistent with the experimental structures of single guide RNAs (sgRNAs). Other clusters, for which there are not yet structures available, were predicted to form novel tracrRNA folds. These clusters included a large and divergent tracrRNA set from Bacteroidetes. These computational models contribute to the understanding of CRISPR-Cas biology, and will assist in the design of further engineered CRISPR-Cas9 systems. The tracrRNA prediction software is available through a galaxy web server.

'Stop' in protein synthesis is modulated with exquisite subtlety by an extended RNA translation signal.

Translational stop codons, UAA, UAG, and UGA, form an integral part of the universal genetic code. They are of significant interest today for their underlying fundamental role in terminating protein synthesis, but also for their potential utilisation for programmed alternative translation events. In diverse organisms, UAA has wide usage, but it is puzzling that the high fidelity UAG is selected against and yet UGA, vulnerable to suppression, is widely used, particularly in those archaeal and bacterial genomes with a high GC content. In canonical protein synthesis, stop codons are interpreted by protein release factors that structurally and functionally mimic decoding tRNAs and occupy the decoding site on the ribosome. The release factors make close contact with the decoding complex through multiple interactions. Correct interactions cause conformational changes resulting in new and enhanced contacts with the ribosome, particularly between specific bases in the mRNA and rRNA. The base following the stop codon (fourth or +4 base) may strongly influence decoding efficiency, facilitating alternative non-canonical events like frameshifting or selenocysteine incorporation. The fourth base is drawn into the decoding site with a compacted stop codon in the eukaryotic termination complex. Surprisingly, mRNA sequences upstream and downstream of this core tetranucleotide signal have a significant influence on the strength of the signal. Since nine bases downstream of the stop codon are within the mRNA channel, their interactions with rRNA, and r-proteins may affect efficiency. With this understanding, it is now possible to design stop signals of desired strength for specific applied purposes.

The exon-intron gene structure upstream of the initiation codon predicts translation efficiency.

Introns in mRNA leaders are common in complex eukaryotes, but often overlooked. These introns are spliced out before translation, leaving exon-exon junctions in the mRNA leaders (leader EEJs). Our multi-omic approach shows that the number of leader EEJs inversely correlates with the main protein translation, as does the number of upstream open reading frames (uORFs). Across the five species studied, the lowest levels of translation were observed for mRNAs with both leader EEJs and uORFs (29%). This class of mRNAs also have ribosome footprints on uORFs, with strong triplet periodicity indicating uORF translation. Furthermore, the positions of both leader EEJ and uORF are conserved between human and mouse. Thus, the uORF, in combination with leader EEJ predicts lower expression for nearly one-third of eukaryotic proteins.

Know Your Enemy: Successful Bioinformatic Approaches to Predict Functional RNA Structures in Viral RNAs.

Structured RNA elements may control virus replication, transcription and translation, and their distinct features are being exploited by novel antiviral strategies. Viral RNA elements continue to be discovered using combinations of experimental and computational analyses. However, the wealth of sequence data, notably from deep viral RNA sequencing, viromes, and metagenomes, necessitates computational approaches being used as an essential discovery tool. In this review, we describe practical approaches being used to discover functional RNA elements in viral genomes. In addition to success stories in new and emerging viruses, these approaches have revealed some surprising new features of well-studied viruses e.g., human immunodeficiency virus, hepatitis C virus, influenza, and dengue viruses. Some notable discoveries were facilitated by new comparative analyses of diverse viral genome alignments. Importantly, comparative approaches for finding RNA elements embedded in coding and non-coding regions differ. With the exponential growth of computer power we have progressed from stem-loop prediction on single sequences to cutting edge 3D prediction, and from command line to user friendly web interfaces. Despite these advances, many powerful, user friendly prediction tools and resources are underutilized by the virology community.

Interference-driven spacer acquisition is dominant over naive and primed adaptation in a native CRISPR-Cas system.

CRISPR-Cas systems provide bacteria with adaptive immunity against foreign nucleic acids by acquiring short, invader-derived sequences called spacers. Here, we use high-throughput sequencing to analyse millions of spacer acquisition events in wild-type populations of Pectobacterium atrosepticum. Plasmids not previously encountered, or plasmids that had escaped CRISPR-Cas targeting via point mutation, are used to provoke naive or primed spacer acquisition, respectively. The origin, location and order of spacer acquisition show that spacer selection through priming initiates near the site of CRISPR-Cas recognition (the protospacer), but on the displaced strand, and is consistent with 3'-5' translocation of the Cas1:Cas2-3 acquisition machinery. Newly acquired spacers determine the location and strand specificity of subsequent spacers and demonstrate that interference-driven spacer acquisition ('targeted acquisition') is a major contributor to adaptation in type I-F CRISPR-Cas systems. Finally, we show that acquisition of self-targeting spacers is occurring at a constant rate in wild-type cells and can be triggered by foreign DNA with similarity to the bacterial chromosome.

Infiltration-RNAseq: transcriptome profiling of Agrobacterium-mediated infiltration of transcription factors to discover gene function and expression networks in plants.

BACKGROUND: Transcription factors (TFs) coordinate precise gene expression patterns that give rise to distinct phenotypic outputs. The identification of genes and transcriptional networks regulated by a TF often requires stable transformation and expression changes in plant cells. However, the production of stable transformants can be slow and laborious with no guarantee of success. Furthermore, transgenic plants overexpressing a TF of interest can present pleiotropic phenotypes and/or result in a high number of indirect gene expression changes. Therefore, fast, efficient, high-throughput methods for assaying TF function are needed. RESULTS: Agroinfiltration is a simple plant biology method that allows transient gene expression. It is a rapid and powerful tool for the functional characterisation of TF genes in planta. High throughput RNA sequencing is now a widely used method for analysing gene expression profiles (transcriptomes). By coupling TF agroinfiltration with RNA sequencing (named here as Infiltration-RNAseq), gene expression networks and gene function can be identified within a few weeks rather than many months. As a proof of concept, we agroinfiltrated Medicago truncatula leaves with M. truncatula LEGUME ANTHOCYANIN PRODUCITION 1 (MtLAP1), a MYB transcription factor involved in the regulation of the anthocyanin pathway, and assessed the resulting transcriptome. Leaves infiltrated with MtLAP1 turned red indicating the production of anthocyanin pigment. Consistent with this, genes encoding enzymes in the anthocyanin biosynthetic pathway, and known transcriptional activators and repressors of the anthocyanin biosynthetic pathway, were upregulated. A novel observation was the induction of a R3-MYB transcriptional repressor that likely provides transcriptional feedback inhibition to prevent the deleterious effects of excess anthocyanins on photosynthesis. CONCLUSIONS: Infiltration-RNAseq is a fast and convenient method for profiling TF-mediated gene expression changes. We utilised this method to identify TF-mediated transcriptional changes and TF target genes in M. truncatula and Nicotiana benthamiana. This included the identification of target genes of a TF not normally expressed in leaves, and targets of TFs from other plant species. Infiltration-RNAseq can be easily adapted to other plant species where agroinfiltration protocols have been optimised. The ability to identify downstream genes, including positive and negative transcriptional regulators, will result in a greater understanding of TF function.

CRISPRDetect: A flexible algorithm to define CRISPR arrays.

BACKGROUND: CRISPR (clustered regularly interspaced short palindromic repeats) RNAs provide the specificity for noncoding RNA-guided adaptive immune defence systems in prokaryotes. CRISPR arrays consist of repeat sequences separated by specific spacer sequences. CRISPR arrays have previously been identified in a large proportion of prokaryotic genomes. However, currently available detection algorithms do not utilise recently discovered features regarding CRISPR loci. RESULTS: We have developed a new approach to automatically detect, predict and interactively refine CRISPR arrays. It is available as a web program and command line from bioanalysis.otago.ac.nz/CRISPRDetect. CRISPRDetect discovers putative arrays, extends the array by detecting additional variant repeats, corrects the direction of arrays, refines the repeat/spacer boundaries, and annotates different types of sequence variations (e.g. insertion/deletion) in near identical repeats. Due to these features, CRISPRDetect has significant advantages when compared to existing identification tools. As well as further support for small medium and large repeats, CRISPRDetect identified a class of arrays with 'extra-large' repeats in bacteria (repeats 44-50 nt). The CRISPRDetect output is integrated with other analysis tools. Notably, the predicted spacers can be directly utilised by CRISPRTarget to predict targets. CONCLUSION: CRISPRDetect enables more accurate detection of arrays and spacers and its gff output is suitable for inclusion in genome annotation pipelines and visualisation. It has been used to analyse all complete bacterial and archaeal reference genomes.